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I.  Theoretical  Phase 
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The  first  summer  of  the  contract  period  was  spent  at  NBS 
Boulder.  Work  done  in  collaboration  with  R.  Peterson  of  NBS 
Boulder  was  published  some  time  ago.  (Phys.  .Rev.  B18,  1198 
(1978).  A  copy  is  attached.  f  j_  f\/  (L  :t'  '■)*■/-<•  '  £  7 

A  successful  program  for  applying  the  microscopic  theory 
to  the  computation  of  the  time-dependent  behavior  of  a  Josephson 
junction  has  now  been  developed.  It  was  based,  originally, 
on  one  due  to  D.  MacDonald  of  NBS  Boulder,  but  extensive  modifica¬ 
tions  had  to  be  made  in  order  to  bring  the  computation  time 
and  computer  memory  requirements  down  to  levels  that  made  it 
possible  to  deal  with  the  kinds  of  problems  I  wished  to  attack. 

Qc  <-  «  ■  t  re*  •  u . 

Such  computations  are  very  cfcrnTandrnq  because  a  iunction- 
has  a  long  "memory".  Its  behavior  at  some  moment  is  dependent 
on  its  history  over  a  long  time  interval  in  the  past.  It  can 
be  shown  that  the  length  of  this  metacry  is  directly  i elated 
to  the  sharpness  of  the  voltage  jump  that  appears  in  tin  junction 
current-voltage  character ist  ic  at  the  gap  voltage .  R.  Harris 
of  NBS  Boulder  has  pointed  oat  that  real  junctions  do  not  have 
as  sharp  a  jump  as  theory  predicts,  and  this  implies  that  their 
memories  are  not  as  long  as  theory  suggests.  1  investigated 
this  question  and  found  that  artificially  shortening  the  theo¬ 
retical  memory  does  indeed  lead  to  rounding  of  the  jump  at 
the  gap  voltage.  1  found  that  the  amount  of  shortening  that 
is  required  to  give  a  reasonable  rounding  is  in  fact  enough 
to  produce  an  appreciable  reduction  in  both  computation  time 
and  computer  memory  requirements . ^  Thi s  change  and  various 
improvements  in  the  program  allowed  accurate  microscopic  theory  ■ 
computati  ns  to  be  made  for  times  as  long  as  160T,  where  T 
is  the  characteristic  gap  time  given  by  Planck's  constant  divided 
by  the  gap  energy  2 A.  I  believe  that  no  one  has  done  anything 
close  to  this  up  to  now. 

At  the  end  of  this  report  I  have  appended  a  copy  of  the 
article  (IEEE  Trans.  Magnetics  MAG- 1 7 ,  809  (1981)  that  describes 
my  calculations,  using  the  above  program,  of  flux  entry  into 
a  single  junction  SQUID.  This  was  the  problem  that  this  contract 
was  intended  to  study.  The  results  bear  out  the  prediction, 
made  in  the  proposal  for  this  project,  that  at  high  frequencies 
the  simple  resistively-shunted  junction  model  is  not  satisfactory 
for  calculating  this  flux  entry. 
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Also  appended  is  another  article  (J.  Appl .  Physics,  accepted 
for  publication)  on  the  remainder  of  the  work  that  I  have  done 
with  the  new  program.  It  gives  some  of  the  details  of  the  method 


the  artificial  "forgetting"  that  I  mentioned  earlier,  and  a  study 
of  the  details  of  the  switching  of  a  junction  from  the  zero- 
voltage  to  non-zero  voltage  state.  As  with  the  SQUID  computations, 
I  found  that  the  RSJ  model  greatly  overestimates  the  amount  of 
damping  that  occurs. 


Experimental  Phase 


1! 
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I  was  forced  to  abandon  the  experimental  part  of  the  project. 
All  of  the  students  who  were  involved  have  loft  Buffalo.  In 
particular,  one  student,  a  Ph.D.  candidate  who  had  worked  for 
two  years,  dropped  out  of  school  for  personal  reasons. 

III.  Applications 

The  high-speed  Joseph  son  junction  circuits  th.it  have  been 
investigated  with  the  program  developed  in  this  project  are 
not  very  different  from  some  proposed  real  device.,,  and  the 
quantities  that  have  been  calculated  are  oi  considerable  interest 
for  device  design.  The  basic  finding  is  that  the  " resist ively 
shunted  junction"  model  that  is  commonly  used  for  device  analysis 
will  in  some  cases  give  poor  results.  Even  fairly  sophisticated 
improved  versions  of  that  model  will  bo  poor  in  some  circumstances 

Further  applications  of  the  new  program  are  now  needed 
in  order  to  determine  when  it  is  important  to  use  it  for  device 
design  and  when  the  simpler  models  previously  used  are  adequate. 
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Smith  and  Blackburn  and  Blackburn  at  al  measured 


flux  entry  Into  single  junction  SOUIDs  and  found 
results  that  disagreed  with  their  calculations,  which 
were  based  on  the  reslstlvely  shunted  junction  (RSJ) 
model.  Since  this  model  Is  known  to  be  very  crude,  it 
seems  worthwhile  to  see  if  use  of  a  more  accurate 
theory  could  remove  the  discrepancy.  Also,  this 
system,  or  a  more  complex  version  of  it,  holds  promise 
for  use  as  a  high  speed  computer  element,  so  It  Is 
Important  to  understand  Its  behavior  In  some  detail. 
This  paper  presents  the  first  study  of  this  SQUID  using 
the  full  BCS  theory  of  superconductivity.  The  results 
are  quite  different  from  the  RSJ  model  predictions,  but 
che  disagreement  with  experiment  Is  made  worse. 


The  experiments  dealt  with  a  large  superconducting 
loop  of  inductance  L  Interrupted  by  a  Josephson  tunnel 
junction.  A  magnetic  field,  applied  perpendicular  to 
the  plane  of  the  loop,  was  slowly  Increased  until  a 
critical  point  was  reached  at  which  a  large  amount  of 
magnetic  flux  would  abruptly  break  Into  the  loop.  The 

measured  quantity  was  6  ,  which  we  will  denote  6  , 

’  enter  a 

the  flux  In  the  loop  at  the  conclusion  of  this  break-in. 


The  equivalent  circuit  for  the  specimen  would  be 
an  inductance  L  In  parallel  with  a  tunnel  junction  of 
critical  current  I  ,  energy  gap  2A,  and  normal  state 


resistance  and  also  In  parallel  with  the  junction 

capacitance  C.  All  wires  are  taken  as  superconducting 
so  that  the  only  losses  would  be  from  the  tunneling 
process  and  so  that  the  magnetic  fluxold  In  the  loop 
would  be  quantized  In  any  steady  state  situation.  In 
the  RSJ  model  the  tunneling  Is  represented  by  a  element 
obeying  the  Josephson  relation  I  «  I  slnS,  where  I  is 

the  supercurrent ,  and  0  Is  the  quantum  mechanical  phase 
difference  across  the  Junction,  In  parallel  with  a  con¬ 
stant  resistance  which  represents  dissipative  quasi- 

particle  tunneling.  In  the  present  work,  on  the  other 

hand,  the  tunnellna  Is  described  using  Werthamer's^ 
analysis  of  the  t lrae-d-. ,  undent  behavior  of  a  Junction 
according  to  the  BCS  theory.  Within  the  contrxt  of 
that  theory,  this  Is  a  completely  general  treatment  of 
the  problem  even  when  the  voltage  la  tlre-dependent. 


I  have  studied  this  system  using  a  computer  simu¬ 
lation  which  is  based  on  the  following  time  domain 

4 

equations  developed  by  Harris. 


where 


I(t)  -  ill  (U*(t)  j”  I^p(t-t")U(t')dt' 
♦  U(t)  j“  Ij(t-t')  U(t')dt'), 

-  exp(-l^  [V,-.*  -  |  1  8  J. 


«> 


U(t) 


(2) 


Here  I(t)  Is  the  tunneling  current  at  time  t  In  terms 


of  the  Junction  voltage  V(t).  The  constant  0g  is  the 


t  "  0  value  of  the  phase  difference  across  the  junction, 
and  ft  Is  assumed  that  V  la  zero  for  t  <  0.  The  ker¬ 
nels  I  (t)  and  I,(t)  are  zero  for  negative  argument, 
qp  J 

For  Identical  superconductors  at  absolute  zero  and  for 
t  »  0  they  are 


I  (t) 
W 


a  ad 


Ij(t) 


*  -2.(fi/eRN)d'(t) 

(3) 

jo  o  yo  <?>• 

(*) 

where  J  and  Y  are  the  nth-order  Bessel  functions  of 
n  a 

the  first  and  second  kind,  respectively.  The  details 
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of  the  method  are  given  elsewhere,  but  one  unusual 
feature  needs  to  be  explained  again  here.  1  have 
Introduced  an  artificial  "memory  reduction"  Into  the 
kernels  I  (t)  and  I  (t)  by  multiplying  each  Bessel 
RP  2  J  2 

function  by  exp(-t  /2At  ),  where  t  •  h/2.*  and  A  Is  a 
constant  which  for  tin  Junctions  Is  taken  to  be  50.  In 
addition,  the  Bessel  functions  were  set  equal  to  zero 
whenever  the  t  In  their  arguments  exceeded  20t .  Ref¬ 
erence  5  shows  that  these  modifications  chance  the 
results  by  no  more  than  a  few  percent  (and  the  change 
■ay  actually  be  an  Improvement  for  real  junctions). 
Their  effect  Is  to  artificially  intioduce  a  small 
amount  of  energy  gap  rounding,  and  the  constant  A  was 
chosen  to  give  the  amount  of  rounding  appropriate  for 
tin.  These  changes  greatly  reduce  the  demands  on  the 
computer  whenever  It  is  necessary  to  follow  the  junc- 
tlonk  behavior  for  a  long  period  ot  tlDe,  and  this 
makes  calculations  feasible  which  otherwise  would  he 
Impractical. 

The  BCS  relation  between  Ic,  R^  and  A  was  assumed 

to  hold,  so  the  only  parameters  that  can  vary  from 
SQUID  to  SQUID  are  L,  I£ ,  R^  and  C.  They  will  be 

expressed  In  terms  of  the  dimensionless  parameters 
t  5  rtc/RjjC  and  y  z  where  Is  the  flux  quan¬ 

tum.  B  is  the  familiar  damping  constant  of  the  LCR 
oscillator,  and  y  is  the  number  of  flux  quanta  pro¬ 
duced  by  a  current  I  circulating  In  the  loop.  The 


basic  Junction  parameter  used  In  reference  5,  R^C/t, 
can  be  found  from  R^C / t  *  2y/«S*. 


One  result  of  the  computations  Is  the  dependence 


of  4  on  y  and  6.  This  Is  shown  in  Fig.  1,  where 


the  dimensionless  quantity  8^/y#o  1*  plotted  against  8. 


Only  cases  for  which  t  /y*  was  less  than  unity  are 

*  C  0 


shown.  For  larger  values,  the  tesults  will  tend  to  be 

6 

erratic  and  very  sensitive  to  initial  conditions. 

The  figure  shows  two  RSJ  points.  They  are  for  y  *  5*0 
but  more  closely  resemble  the  BCS  reiults  for  y  •  '0. 
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Figure  1.  Dependence  of  the  amount  of  flux  entering 


the  SQUID  on  the  damping  parameter  8*  y*o  la  LIc. 

All  of  the  points  except  the  solid  circles  were  com¬ 
puted  using  the  BCS  theory.  The  solid  circles  wera 
taken  from  Fig.  2  of  reference  1. 

To  shift  one  of  these  points  to  the  BCS  value  would 
require  Increasing  the  6  value  used  In  the  graph, 
which  could  be  regarded  as  decreasing  the  effective 
resistance  urcd  In  calculating  B.  On  the  other  hand. 
If  one  calculates  0  from  R^,  then  the  correct  could 

be  obtained  by  using  a  larger  "effective  R^"  in  the 

computer  simulation,  which  would  reduce  the  amount  of 
damping  and  give  a  larger  6^.  In  this  sense,  the 

original  RSJ  calculation  overestimates  the  amount  of 
damping. 

In  addition,  some  qualitative  conclusions  from  the 
RSJ  analysis  do  not  apply  in  the  BCS  case.  Smith  and 

Blackburn1  found  an  expression  for  a  quantity  Pc(y) 
which  gave  the  minimum  6  such  that  6^  was  one  quantum. 
This  same  0£  gave  a  universal  curve  when  6^/y  was 
plotted  against  0/0£.  There  seems  to  be  no  equivalent 

quantity  for  the  BCS  results.  In  fact,  for  at  least 
one  ease  (y-50)  there  6eems  to  be  no  sensible  value  of 
0  that  will  restrict  to  one  quantum.  In  Fig.  I, 

4  /yt  seems  to  be  approaching  a  limit  of  about  0.1 
*  o 

(*  -  56  ).  Even  with  the  totally  unrealistic  value 

0  O 

of  B  of  5600  (not  shown)  ♦  It  still  5#  . 

s  o 

An  Interesting  feature  was  observed  for  thl#  caae 
of  y-50, 8  •  5600.  The  time  derivative  of  the  flux  In 
the  loop  became  slightly  negative  briefly  when  wa» 

1.25«o.  It  then  became  positive  again,  and  con¬ 
tinued  to  increase  to  5.36  .  This  It  quite  inexplic¬ 
it 

able  by  the  RSJ  model  (even  If  modified  to  let  the 
resistance  be  voltage  dependent).  In  that  picture, 
once  the  flux  hat  started  to  decrease  In  the  potential 
energy  well  corresponding  to  a  particular  number  of 
quanta  the  system  will  be  trapped  In  that  well. 

Blackburn  et  al2  proposed  that  the  principle 


source  of  disagreement  between  their  measurements  and 
their  RSJ  calculation  was  the  effects  of  finite  quasl- 
partlcle  lifetimes  In  the  metal  films  of  the  Junctions. 
1  have  not  actually  calculated  the  behavior  of  Smith 

and  Blackbutn'a  specimen,  which  had  y  •  5  x  10*.  Even 
with  memory  reduction,  such  a  large  value  of  y  makea 
the  computer  simulation  impractical.  However,  because 
of  the  univcrssl  curve  for  6^/y  found  by  them,  there 

Is  no  reason  to  doubt  that  for  this  case  the  BCS  and 
RSJ  calculations  will  differ  In  the  same  general  way  as 
they  do  for  smaller  y.  In  other  words.  It  Is  clear 
that  the  BCS  calculation  will  give  greater  disagreement 
with  the  experiment,  and  it  seems  unlikely  that  finite 
lifetime  effects  can  be  responsible. 

For  device  considerations  It  it  also  of  interest 
to  determine  the  time  required  for  flux  entry  in  the 
SQ’JID.  There  is  an  Initial  period  of  very  slow  change, 
until  the  flux  Is  near  one  quantum.  For  the  rest  of 
the  tine  the  average  junction  voltage  is  close  to  the 
gap  voltage,  and  flux  enters  at  the  corresponding  rate 

of  26  ft.  The  time  for  the  first  quantum  to  enter, 

o 

which  I  will  call  t^. depends  only  on  R^C  for  the  range 
of  parameters  studied.  As  Fig.  2  shows,  t^/r 
approaches  a  minimum  value  of  4.7  for  small  R^C.  For 
larger  R„C,  It  becomes  roughly  proportional  to 


(RjjC/t)5. 


(RnC/  T 


Figure  2.  Dependence  on  device  parametera  of  the  time 
t^  for  the  first  quantum  of  flux  to  enter  the  SQUID, 

according  to  the  BCS  theory.  When  t./r  la  plotted 

w  1 

againat  (R^C/t)  ,  a  universal  curve  results  that  la 
roughly  linear  for  larger  (The  8  values  for  each 

point  esn  be  computed  from  8  “  (2yt/  xV)Y 

The  RSJ  model  gives  similar  behavior,  except  that 
the  average  voltage  falla  appreciably  below  the  gap 
value  in  the  later  stages  of  flux  entry.  As  an  example, 
our  own  unpublished  RSJ  calculations  give  the  following 
results  for  y  •  120,  8  "  5:  is  120r,  and  the  total 

time  to  reach  the  final  value  of  *  936^  la  238i,  ao 


f 


V  ' 


V 


■» 

that  d+/dt  averaged  .8$  h  after  t,.  The 
O  1 

BCS  value*  are:  t,  ■  26r,  ♦  *  1244  ,  total  time 
1  8  0 

-  88t,  and  the  average  rate  after  t,  la  2.04  /t.  Sailth 

1  o 

and  Blackburn  did  not  measure  theae  tinea  (they 
are,  of  course,  extremely  short),  so  It  is  not  pos- 
sible  to  compare  theory  and  experiment.  In  any 
case,  ve  again  aee  that  the  RSJ  nodel  overeat laates 
the  damping. 

In  summary,  the  simple  constant  resistance  RSJ 
model  gives  poor  results  for  flux  entry  into  a 
alngle-Junct ton  SQUID  and  is  not  suitable  for  quanti¬ 
tative  comparison  with  experiment.  (A  modified  RSJ 
model  in  which  a  kind  of  energy  gap  is  Introduced 
by  giving  the  resistance  a  suitable  voltage  de¬ 
pendence  may  work  better,  but  computer  simulation 

of  single  Junction  switching  shows  that  even  this 
is  not  necessarily  satisfactory.)  The  predictions 
of  the  BCS  theory  do  not  agree  with  experiment. 

The  reason  for  this  is  not  known. 
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Josephson  tunnel  junction  time-dependence 
calculations  based  on  the  BCS  theory 
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ABSTRACT 

The  feasibility  and  the  importance  of  using  the  full  BCS 
theory  for  numerical  calculations  of  the  time-dependent  behavior 
of  a  Josephson  tunnel  junction  is  examined  by  investigating 
switching  from  the  zero  to  the  finite  voltage  state.  The  cal¬ 
culation  can  be  done,  but  only  if  the  RC  time  constant  is  not 
too  long.  The  results  differ  appreciably  from  those  obtained 
with  the  crude  but  frequently  used  "resistively  shunted  junction 
model.  Artificial  rounding  of  the  energy  gap  edge,  in  the  form 
of  "memory  reduction",  is  introduced  and  shown  to  reduce 
computational  difficulties  appreciably  without  significantly 
altering  any  results.  This  should  make  many  more  situations 
amenable  to  accurate  digital  simulation. 
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I.  Introduction 

The  qualitative  behavior  of  a  Josephson  tunnel  junction 
can  be  understood  using  the  very  simple  "resistively  shunted 
junction"  (RSJ)  model,  but  accurate  calculations,  particularly 
at  high  frequencies  or  where  damping  is  important,  require  a 
more  detailed  analysis.  This  paper  presents  the  start  of  an 
investigation  of  the  consequences  of  applying  the  BCS  micro¬ 
scopic  theory  to  time-dependent  problems.  We  will  see  that 
in  some  instances  the  quantitative  results  are  quite  different 
from  those  obtained  with  the  RSJ  model. 

The  analysis  has  been  done  by  computer  simulation,  and 
the  method  used  will  be  discussed  in  the  following  section. 

In  addition,  some  computations  using  a  modified  form  of  the 
theory  which  gives  nearly  the  same  results  but  greatly  reduced 
computation  cost  will  be  described. 

II.  Method  of  Calculation 

The  BCS  theory  was  used  to  calculate  junction  behavior 

1  2 
by  applying  Harris'  time  domain  formulation  of  Werthamer's 

analysis  of  junction  dynamics.  •  The  basic  relation  gives  the 

tunneling  current  I  at  time  t  in  terms  of  the  voltage  V  at  all 

times  up  to  and  including  t.  In  Harris'  notation  it  is 


I(t) 


Im  (U*(t) 


°o 

I  (t-t')u(t ')dt' 
qp 

—  oo 


+U(t) 

where 


I j (t-t ' ) U (t '  )dt '  1 , 

U(t)  =  exp  l-i  r 


V(t')dt'  +  ~  i  <|>0] 


(la) 


(lb) 
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The  constant  <t>Q  is  the  t  =  0  value  of  the  phase  difference  across 

the  junction,  and  it  is  assumed  that  V  is  zero  for  t  <  0. 

Equation  (la)  appears  to  require  knowledge  of  V  for  all  times 

in  the  future  as  well  as  the  past,  but  this  is  not  so  since  the 

kernels  I  (t)  and  I  (t)  are  zero  for  negative  argument.  For 
^IP  J 

identical  superconductors  at  absolute  zero  and  for  t  >  0  they  are 


I  <t)  =  —  2  tt  (h/eRN)«'(t) 

V«  -js¥Jo  (T>,o‘lf>-  <2» 

N 

Here  is  the  normal  state  tunneling  resistance  of  the  junction, 

2A  is  the  energy  gap,  and  and  are  the  nth-order  Bessel 
functions  of  the  first  and  second  kind,  respectively. 

This  represents  an  exact  microscopic  description,  based 
on  the  BCS  theory,  with  the  following  limitations:  (1)  Only 
identical  superconductors  at  absolute  zero  are  considered  here; 

(2)  The  junction  must  be  small  so  that  spatial  variations  can 
be  ignored;  (3)  No  nonequilibrium  or  quasiparticle  lifetime 
effects  are  included;  (4)  Energy  gap  anisotropy  is 
ignored;  (5)  No  losses  in  the  superconducting  metal  are  included. 

A  second  equation  relating  1 { t )  and  V ( t )  will  be  given  by 
the  specific  circuit  being  considered.  If  the  initial  values  <J> 
and  V(0)  are  given,  the  two  equations  can  be  integrated  numerically 
to  give  I (t)  and  V(t).  (The  junction  capacitance  C  is  treated 
as  part  of  the  external  circuit.) 


«Wl mm''** ' 


Harris  used  this  approach  to  find  the  effect  of  applying  a 
voltage  pulse, ^  to  examine  the  switching  of  a  junction  having  an 
external  load  resistance,  and  to  investigate  the  small,  sinusoidal 
voltage  limit.  ^  McDonald,  Johnson,  and  Harris'*  analyzed  a  few 
current-biased  cases  to  check  their  frequency  domain  calculations 
of  current-voltage  characteristics.  Noting  that  Eq .  (la)  requires 
V  at  time  t  to  determine  I  at  the  same  time,  and  that  (la)  is 
nonlinear  in  V  because  of  the  form  of  U,  they  used  an  iteration 
scheme  that  starts  with  a  guessed  value  for  V.  In  the  present 
calculations  a  simpler  procedure  has  been  used.  Equation  (lb) 
was  expanded  in  a  power  series  to  first  order  in  the  quantity 

(At/x) (V  -V  i , ) /V  . ,  ,  where  x  is  h/2A.  This  is  of  order 
new  old  old 

2 

(At/x )  ,  which  is  the  order  to  which  all  approximations  were 

made.  Equation  (la)  is  thereby  linearized  and  can  be  solved 
simultaneously  with  the  linear  circuit  equation. 

The  accuracy  of  the  results  was  investigated  by  comparing 
computations  done  with  and  without  double  precision  arithmetic, 
by  varying  the  time  step  size,  and,  where  possible,  by 
comparing  with  results  of  frequency-domain  calculations. 

These  tests  indicated  that  all  results  presented  here  are  good 
to  2  percent  or  better. 

Anyone  interested  in  further  details  of  the  computation 
method  is  urged  to  contact  the  author  directly.  However,  it  is 
important  to  note  that  the  form  of  Eqs.  (1)  and  (2)  make  the 
calculation  very  demanding.  First  of  all,  the  kernels  in 
Eq.  (1)  decay  as  t  \  which  is  very  slow.  The  integrals  in 
Eq.  (la)  are  of  the  form  of  convolutions  and  must  be  done 
all  over  again  for  each  time  step.  This  means  that  the  com¬ 
putation  time  will  increase  as  t2.  it  also  means  that 


every 
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value  of  I  ,  I_,  and  the  real  imaginary  parts  of  U  must  be 
qp  J 

remembered.  Depending  on  the  computing  facilities  available, 
either  the  computing  time  or  the  memory  requirement  can  limit 
the  maximum  t  that  can  be  considered. 

These  difficulties  are  compounded  by  two  facts.  One  is 
that  typical  junctions  have  RC  time  constants  of  the  order  of 
10t  or  more.  The  other  is  that  x,  a  characteristic  time  for  a 
junction,  is  of  the  order  of  picoseconds,  which  is  very  short 
compared  to  the  time  scale  of  most  measurements.  The  second 
fact  means  that  we  sometimes  cannot  numerically  duplicate  the 
procedure  of  an  experiment,  such  as  the  determination  of  a 
current-voltage  characteristic,  but  must  start  with  unphysical 
initial  conditions.  The  first  fact  means  that  the  inevitable 
transients  produced  by  this  will  require  many  x ’ s  to  die  away. 

As  an  example,  consider  the  I-V  characteristic  for  R^C  =  lOx 
which  will  be  presented  below.  Each  point  was  obtained  from  a 
separate  calculation.  The  time  step  necessary  to  get  good 
results  was  .02x,  and  it  was  necessary  to  compute  for  80x  in 
order  to  determine  the  average  voltage  to  within  less  than  1% 
of  its  limiting  value.  Thus,  4000  time  steps  were  required, 
as  well  as  the  storage  of  four  4000-element  arrays.  It  is  there¬ 
fore  not  surprising  that  with  the  existing  computer  program 
it  was  not  feasible  to  work  with  larger  values  of  R^C. 

There  are  various  steps  which  can  be  taken  to  reduce  these 
difficulties.  One  is  to  improve  the  approximations  made  in 
the  numerical  method  and  thereby  allow  longer  time  steps  to 
be  used.  This  is  now  being  worked  on. 


Tier 


Another  modification  is  based  on  the  fact  that  the  kernels 


1^  and  I  do  eventually  become  negligibly  small.  For  large 

enough  t-t '  they  can  be  set  equal  to  zero.  Beyond  this  point  the 

2 

computing  time  goes  as  t  instead  of  t  ,  and  the  size  of  the 
arrays  that  must  be  remembered  does  not  increase.  In  practice 
this  has  not  been  found  helpful  except  when  combined  with 
another  change,  a  reduction  in  the  length  of  the  junction's 


"memory'' . 

"Memory"  refers  to  the  fact  that  I  at  time  t  depends 

on  V  at  all  previous  times.  The  slow  fall  off  of  the  kernels 

3 

gives  a  very  long  memory.  Harris  has  pointed  out  that  the 


length  of  the  memory  is  directly  related  to  the  sharpness 
of  the  rise  in  current  seen  in  the  I-V  characteristic  at  the 


gap  voltage.  The  rise  exhibited  by  real  junctions  is  never  as 
sharp  as  is  predicted  by  the  BCS  theory.  Harris  suggested 
that  this  means  that  the  memory  of  real  junctions  is  shorter 
than  predicted,  or  in  other  words,  the  real  kernels  fall  off 
more  rapidly  than  as  t  \  Perhaps  a  more  detailed  theory, 
including  such  complications  as  gap  anisotropy,  would  give 
modified  kernels  leading  to  a  shorter  memory  and  thereby 
simplifying  numerical  computations. 

In  lieu  of  such  a  more  detailed  theory,  the  consequences 
of  a  phenomenological  reduction  in  memory  have  been  investigated, 
and  the  results  are  rather  encouraging.  In  some  of  the  cal¬ 
culations  to  be  described  below,  each  of  the  Bessel  functions 
in  Eq.  (2)  was  multiplied  by  exp  (-t  /2At  ).  The  constant  A 
was  chosen  to  be  50,  which  altered  the  I-V  characteristic  by  a  few 
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percent,  giving  a  rounding  that  matched  fairly  closely 
what  was  observed  in  tin  tunnel  junctions.^  Instantaneous 
current  and  voltage  values  were  also  changed  by  only  a 
few  percent,  and  usually  less  than  two  percent,  which 
for  most  practical  purposes  is  not  significant.  This 
allows  one  to  hope  that  a  modification  based  on  physical 
grounds  would  also  not  alter  the  detailed  time  dependence 
by  a  large  amount  and  that  the  Werthamer  theory  as  well 
as  the  phenomenological  modification  introduced  here  will 
not  disagree  substantially  with  a  more  complete  treatment. 

In  any  case,  the  memory  reduction  used  does  not  degrade 
the  results  significantly,  and  it  has  the  great  advantage 
that  it  allows  the  kernels  to  be  replaced  by  their  limiting 
value  of  zero  at  much  smaller  t.  As  explained  earlier, 
this  greatly  reduces  the  demands  made  on  the  computer. 

For  R^C  =  10 T  and  A  =  50,  this  replacement  could  be  made 
at  t  *  20  t. 

Besides  giving  a  more  realistic  I-V  characteristic,  this 
memory  reduction  will  be  shown  to  reduce  and  broaden  the  Riedel 
peak,^  which  is  what  one  would  expect.  A  side  effect  is  a 
reduction  in  Ic,  the  maximum  zero-voltage  supercurrent.  It 
is  not  known  whether  this  is  physically  correct  or  net,  but 
since  the  change  is  only. 0042!  for  the  A  value  used,  it  is 
of  no  practical  consequence. 
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One  goal  of  this  work  is  to  examine  the  accuracy  of  the 
RSJ  model  by  comparing  its  predictions  with  those  of  the  BCS 
theory.  The  basic  RSJ  model  treats  the  junction  as  a  resistor 
R,  equal  to  the  normal  state  resistance  R^,  in  parallel  with 
a  capacitor  C  and  also  in  parallel  with  an  element  obeying 
the  equation  I  =  Icsim-,  where  is  the  superconducting  phase 
difference  across  the  junction.  There  is  nothing  corresponding 
to  an  energy  gap  in  this  picture,  and  so  of  course  the  resulting 
I-V  characteristic  shows  no  structure  at  the  gap. 


m 
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There  is  an  extensive  literature  of  calculations  of  junction 
behavior  with  this  model,  and  many  features  can  be  understood 
qualitatively  in  this  way.  The  simplest  and  probably  most 
important  improvement  that  can  be  made  in  the  model  is  to 
let  R  be  voltage  dependent.  This  can  obviously  give  a  fairly 
realistic  looking  I-V  characteristic  if  R(V)  is  chosen 
appropriately,  and  it  probably  gives  greatly  improved  results 
for  any  calculation.  Results  obtained  with  both  the  basic  RSJ 
model  and  the  variable  resistance  version  will  be  presented. 

In  the  latter,  R(V)  was  chosen  so  that  I  =  V/R  would  fit  our 
experimental  I-V  curves  for  tin  junctions^.  These  junctions 
had  an  Rj^C/t  value  of  about  60,  rather  than  the  10  that  is  used 
in  most  of  these  calculations.  However,  the  I-V  curve  is  not 
expected  to  change  significantly  in  this  particular  range 
of  RnC/t.  From  the  description  of  the  equivalent  circuit, 
one  can  see  that  the  equation  that  replaces  Eq.  (1)  is 


I  =  Ic  sin<J>  +  V/R, 


where 


<t>  - 


2e 


V(t')  dt'  +  <{>  . 


(3a) 


(3b) 


The  junction  capacitance  is  treated  as  part  of  the  external 
circuit,  as  before.  Equations  (3)  together  with  the  appropriate 
circuit  equation  can  be  numerically  integrated  in  a  straight¬ 
forward  way.  The  "memory"  or  dependence  on  past  voltage  is  of  a 
simple  nature  in  this  case  and  does  not  lead  to  the  computational 
difficulties  found  with  the  BCS  theory. 
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All  of  the  results  will  be  expressed  in  normalized  units. 

Current  is  divided  by  the  critical  current  I  =  irA^eR^,  and 

voltage  by  the  gap  voltage  Vg  =  2A/e.  The  unit  of  time  is 

t  =  h/2A.  As  before,  R^  is  the  junction  normal  state  resistance 

and  2 A  is  the  energy  gap  at  absolute  zero.  The  only  adjustable 

parameter  that  can  distinguish  one  junction  from  another  is 

the  dimensionless  RC  time  constant  R^C/t ,  which  is  independent 

of  junction  area.  It  is  proportional  to  the  oxide  barrier's 

dielectric  constant  to  thickness  ratio  and  inversely  proportional 

to  the  maximum  supercurrent  per  unit  area.  In  practice,  for  a 

given  material,  the  former  varies  little  from  junction  to  junction, 

but  the  latter  can  sometimes  vary  over  several  orders  of  magnitude. 

For  tin-tin  oxide-tin  junctions  R^C/r  =  10  would  correspond  to 

2  2 

the  reasonable  supercurrent  density  value  of  about  10  A/cm  . 

This  is  the  value  used  in  this  paper,  but  it  should  be  borne 
in  mind  that  many  aspects  of  a  junction's  behavior  depend  rather 
strongly  on  this  parameter. 

III.  Applications 
A.  Current-voltage  characteristics 
Since  voltage  and  frequency  are  directly  proportional 
for  a  Josephson  junction,  the  natural  approach  to  finding  the 
I-V  curves  from  Werthamer's  analysis  is  to  work  in  the  frequency 
domain.  This  has  been  done  by  McDonald  et.al.4  and  by  Schlup.8 
The  time  domain  formulation  is  used  to  compute  i-v  characteristics 
here  in  order  to  check  that  the  program  is  correct  and  to  in¬ 
vestigate  the  phenomenological  memory  reduction  that  was 
described  in  the  previous  section. 
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The  circuit  in  this  case  is  a  Josephson  junction  in 


parallel  with  a  capacitor  C  and  the  combination  in  series 
with  a  constant  current  source.  The  tunneling  current,  the 
displacement  current,  and  the  voltage  will  vary  rapidly  and 
sometimes  with  large  amplitude.  However,  the  discussion  here 
will  be  limited  to  the  average  voltage,  since  this  is  what 
is  measured. 

With  RjjC/t  =  10  a  time  step  of  .  02x  gave  instantaneous 
current  and  voltage  values  to  about  1  percent  or  better.  The 
initial  value  of  phase,  <pQ,  was  chosen  so  that  the  initial 
tunneling  current  equalled  the  current  1^  from  the  constant 
current  source.  The  initial  voltage  was  usually  taken  to  be 
the  gap  value.  This  is  unphysical  since  the  formalism  assumes 
V  =  0  for  all  t  <  0.  This  led  to  rather  violent  transient 
behavior,  and  so  the  computation  was  carried  out  to  large 
enough  t  so  that  the  average  voltage  no  longer  changed.  (Even 
though  it  required  allowing  the  transients  to  die  out,  this 
scheme  seemed  to  require  less  computation  time  than  would 
others  that  could  be  imagined.)  Then  the  process  was  repeated 
with  a  new  value  of  I^c-  In  order  to  check  that  the  maximum 
zero-voltage  current  had  the  correct  value,  some  computations 
with  I(jc  in  the  neighborhood  of  Ic  were  given  an  initial 
voltage  of  zero.  For  I^c  <  I  the  voltage  remained  zero,  as 
it  should.  For  I^c  <  1  but  initial  voltage  equal  to  the  gap 
value,  the  average  voltage  was  nonzero. 

The  result  is  shown  by  the  circles  in  Figure  1.  The 
shape  of  this  curve  is  consistent  with  the  frequency  domain 

A  O 

results  of  McDonald  et.  al.  and  of  Schlup.  The  accuracy 
of  the  voltage  values  is  believed  to  be  1%.  Note  that  the 
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sharp  current  rise  at  the  gap  voltage  has  a  voltage  width 
of  2%,  so  it  is  not  possible  to  be  positive  that  this  width  is 
real.  However,  the  small  foot  on  the  low  voltage  side  of  the 
current  rise  is  real.  The  dashed  line  is  drawn  from  a  point 
computed  for  Idc/*c  “  10  and  haa  a  slope  corresponding  closely 

to  i^. 

The  solid  line  is  taken  from  I-V  curves  measured  on  our 

high-quality  tin  tunnel  junctions.6  Notice  that  the  current 

rise  at  the  gap  voltage  is  not  as  sharp  as  theory  predicts. 

For  tin,  there  is  evidence  that  this  is  due  to  the  large  variation 

9  10 

in  the  size  of  the  gap  for  different  crystallographic  directions.  ' 
However,  it  seems  that  junctions  of  other,  less  anisotropic, 
materials  show  comparable  rounding.  For  the  present  phenomeno¬ 
logical  approach  it  is  assumed  that  the  physical  cause  of  the 
rounding  is  not  important.  At  worst,  the  memory  red  -tion  introduced 

does  not  nroduce  a  large  change  in  the  results,  and.  at  best  it  may 
actually  improve  them,  just  as  it  improves  the  I-V  curve. 

The  crosses  show  the  change  produced  by  memory  reduction. 

The  parameter  A  was  chosen  to  be  50  by  requiring  the  voltage 
value  at  /I  -  1*0  to  match,  that  observed  in  the  tin  tunnel 
junctions.  The  result  is  very  plausible,  being  closer  to  the 
experimental  curve  than  is  the  unmodified  curve.  It  confirms 
Harris'  statement  that  gap  rounding  will  result  from  memory 
reduction. 
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B.  Constant-voltage  coefficients 
When  the  junction  voltage  is  constant,  Josephson**  showed 
that  the  tunneling  current  can  be  written 

I (V ,T)  =  IJ1(V,T)  sin*  +  Ij2(V,T)  cos *  +  Igp(V,T),(4) 

where  T  is  the  absolute  temperature  and  the  notation  is  that 
12 

used  by  Harris.  In  most  situations  the  voltage  will  not  be 
constant  because  typical  junction  resistances,  of  the  order  of 
an  ohm,  are  much  smaller  than  usual  source  impedances.  On  the 
other  hand,  a  typical  RC  time  constant  of  the  order  of  IOt 
keeps  the  voltage  from  varying  rapidly  in  a  time  interval  t. 


-u- 


Therefore,  there  may  be  situations  in  which  Eq.  (4)  will  give 
fairly  good  results.  It  gives  the  basic  RSJ  model  if  one  sets 
1^^  =  Ic,  1^2  =  0/  and  I  =  V/R^  Obvious  improvements  would 
be  to  give  1^,  Ij2'  an£^  Jqp  t^e  vo*ta9e  dependences  predicted  by 
theory.  Even  then,  however,  the  model  is  strictly  correct  only 
for  constant  voltage.  There  seems  to  be  no  way  to  estimate  its 
precision  except  by  comparison  with  experiment  or  with  the  full 
time-dependent  theory.  In  this  paper  this  process  of  comparison 
will  be  begun  by  examining  the  RSJ  model  and  the  modification  in 
which  I  has  the  voltage  dependence  suggested  by  experiment. 


The  coefficients  of  Eq.  (4)  were  computed  for  T  =  0  using 
the  time-dependent  program  described  earlier.  This  serves  as 
another  check  on  the  program  and  also  allows  determination  of 
the  effect  of  memory  reduction  on  these  coefficients. 

The  circuit  equation  for  this  case  is  simply  the  statement 
that  the  voltage  is  constant.  The  junction  capacitance  p] ays 
no  role.  In  the  program,  the  terms  of  Eq.  (4)  were  not 
identifiable,  so  they  could  not  be  computed  directly.  Instead, 
the  voltage  was  chosen  and  the  computation  was  begun  with  I  and 
4>0  set  to  zero.  The  effect  of  this  is  as  if  a  voltage  step 
function  were  applied  at  t  =  0.  The  program  was  run  until  the 
transients  died  out.  (A  time  step  of  .01t  was  found  to  be 
suitable,  so  the  computation  could  be  carried  out  to  t  =  40t.) 
Then,  from  the  I  and  4>  values  at  three  t’s  the  three  coefficients 
of  Eq.  (4)  were  determined.  This  process  was  repeated  for 
each  voltage. 
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For  the  unmodified  BCS  calculation,  it  was  not  feasible 

to  get  accurate  results  in  the  region  within  ±5%  of  the  gap 

voltage  since  there  the  transients  take  too  long  to  die  out. 

Of  course,  some  such  difficulty  is  inevitable  with  any  method 

because  of  the  singularity  (the  "Riedel  singularity"^)  in  I  ^ 

at  the  gap  voltage.  Outside  of  this  region  the  results  agree 

12 

with  standard  calculations  such  as  those  by  Harris,  whose 
curves  are  reproduced  as  lines  in  Fig.  2. 

The  points  in  Fig.  2  show  the  changes  produced  by  memory 
reduction,  with  A  =  50  as  before.  The  Riedel  peak  is  broadened 
and  reduced,  and  the  abrupt  jumps  in  and  I^p  are 
broadened.  These  are  the  kind  of  alterations  one  would  expect 
from  a  physical  gap-rounding  mechanism,  and  this  lends  support 
to  the  idea  that  the  memory  reduction  introduced  in  this  paper 
may  actually  improve  the  agreement  between  theory  and  experiment, 
not  only  for  I-V  curve  measurement,  but  for  other  cases  as  well. 

C.  Switching 

This  section  deals  with  switching  of  a  junction  from  the 
zero  voltage  to  the  finite  voltage  state  by  a  constant  current 
source.  This  was  chosen  because  it  is  of  practical  interest, 
it  can  be  treated  in  a  fairly  realistic  way,  and  it  is  a  case 
where  the  RSJ  model  gives  rather  bad  quantitative  results. 

The  circuit  is  the  same  as  for  the  I-V  curve  calculation. 

The  initial  conditions  are  V  =  0  and  j;  such  that  1  =  1. 

T  >  C 

For  the  results  presented  here  was  taken  to  be  1.001  Ic. 

These  are  quite  reasonable  initial  conditions.  One  can  imagine 
that  for  the  time  interval  -»  <t<0  the  applied  current  was 
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exactly  Ic  and  the  voltage  remained  zero.  Then  at  t  =  0  a 
small  increase  in  current  caused  switching  to  begin.  This  could 
correspond  to  a  real  experiment  and  should  show  the  true  time 
dependence  of  a  junction  being  used  as  a  switch,  except  for 
certain  approximations  that  have  been  made  in  the  treatment  of 
the  circuit.  The  switching  will  turn  out  to  be  quite  fast, 
so  a  true  constant  current  source  would  be  hard  to  achieve. 

Also,  radiation  will  occur,  but  it  is  not  included  in  the 
equations.  The  latter  is  not  likely  to  be  of  any  consequence 
to  the  switching,  but  the  former  may  well  be,  depending  on  the 
specific  circuit.  The  present  calculation  will  exhibit  the 
behavior  that  is  intrinsic  to  the  junction  itself,  and  there 
would  be  no  particular  difficulty  associated  with  repeating  the 
calculation  with  a  modified  external  circuit. 

The  line  labeled  "BCS"  in  Fig.  3(a)  shows  what  happens 
for  Rj^C  =  10t.  The  other  line  shows  the  P.SJ  result. 

Note  that  while  the  two  calculations  give  qualitatively  similar 
results,  the  details  are  quite. different . 

One  inevitable  failing  of  the  RSJ  model  is  that  it  gives 
the  wrong  final  voltage  because  it  does  not  give  the  correct 
resistance  near  I  =  lc>  This  is  probably  of  no  particular 
concern.  The  switching  time  is  likely  to  be  what  one  needs 
to  predict.  For  this,  the  RSJ  model  fails  in  two  ways.  First, 
it  gives  a  longer  delay  before  the  onset  of  the  rapid  voltage 
rise.  Second,  it  gives  a  more  or  less  exponential  approach 
to  the  final  voltage,  whereas  the  BCS  calculation  shows  a 
steeper,  nearly  linear  rise.  For  example,  the  RSJ  case  rises  from  10% 
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of  its  final  voltage  to  90%  in  a  time  internal  of  20t,  while 
the  BCS  case  rises  from  10%  of  its  final  voltage  to  90%  in 
only  11t.  The  BCS  result  also  differs  from  the  RSJ  in  that 
the  voltage  shows  an  overshoot  at  the  end  of  the  rapid  rise. 

Calculations  using  a  voltage-dependent  resistance  in  the 
RSJ  model  give  much  better  agreement  with  BCS  theory,  as  shown 
in  Fig.  3(b).  For  this  calculation,  the  junction  resistance 
is  very  large  until  the  voltage  gets  very  close  to  the  gap 
value.  Therefore,  during  most  of  the  switching  period  the 
current  is  nearly  all  supercurrent  (lcsin<{>)  or  displacement 
current  through  the  capacitor.  Apparently,  much  the  same  thing 
is  true  for  the  BCS  case.  When  the  rapid  rise  beqins,  the  RSJ  cal 
culation  has  the  voltage  rising  exponentially  toward  the  very 
large  value  that  the  high  resistance  requires.  When  the  voltage 
is  approximately  the  gap  value,  the  resistance  drops  rapidly, 

sharply  cutting  off  the  rapid  voltage  rise.  This  gives  a  fairly 
linear  portion  during  that  rise.  The  overshoot  is  still 

not  reproduced,  however. 

3 

A  similar  comparison  was  carried  out  by  Harris  for  R^C/t  = 
He  considered  the  switching  of  a  junction  that  had  an  external 
load  resistor,  so  the  calculations  cannot  be  compared  directly, 
but  the  qualitative  conclusions  are  the  same. 

The  important  role  played  by  junction  capacitance  can 
be  seen  by  comparing  Fig.  3  with  b,  which  shows  the  R^C  =  0.1  t 

case.  The  switching  is  much  faster  and  the  voltage  oscillations 
are  greatly  enlarged  in  the  latter.  Again,  the  simple  RSJ 
model  (not  shown)  is  in  error.  It  gives  the  first  peak  at 
7.9t  instead  of  3.3t  and  1.2t  instead  of  ,b5T  for  the  oscillation 
period.  More  surprisingly,  the  variable  resistance  case, 
shown  by  dots,  is  also  not  very  successful.  The  explanation 


for  the  disagreement  is  an  interesting  lesson  in  the  hazard  of 
thinking  in  terms  of  Eq.  (4).  As  the  switching  begins,  V 
starts  to  change.  dV/dt  is  very  small  but  not  zero  (and  it  is 
about  100  times  larger  than  for  r^C/t =  10),  and  this  leads 
to  a  small  additional  contribution  to  the  tunneling  current 
besides  what  one  would  expect  from  Eq.  (4) .  In  fact,  in  the 
early  moments  the  tunneling  current  rises  slightly  above  Ic- 
Over  a  period  of  time,  the  cumulative  effect  of  this  additional 
tunneling  current  and  the  corresponding  decrease  in  the 
displacement  current  slows  the  voltage  rise  enough  to  produce 
the  delay  in  switching  compared  to  the  variable  resistance 
RSJ  model  that  Fig.  4  shows. 


These  calculations  were  repeated  with  memory  reduction. 
(The  parameter  A  was  set  to  50,  and  the  kernels  were  cut  off 
for  t  in  their  arguments  of  20t.)  The  results  were  nearly 
identical  to  the  originals,  except  for  a  slight  shift  along 
the  time  axis.  For  R^C  =  10t,  in  the  vicinity  of  the  peak 
at  33t ,  the  shift  was  .34t.  For  R^C  =  O.lx,  around  the  peaks 
near  3t,  the  shift  was  .05r.  It  is  unlikely  that  these 
differences  would  be  of  any  practical  consequence,  so  it  seems 
sensible  to  use  memory  reduction  for  any  further  computations 


of  this  kind. 


In  summary,  one  can  see  that  the  RSJ  model  gives  incorrect 
numbers  and  does  not  really  give  the  right  shape  for  the  voltage¬ 
time  curve.  Simply  letting  R  vary  with  V  can  remove  most  of 
the  discrepancies  for  an  R^C/t  of  10  but  not  for  0.1.  Further 
work  is  needed  to  more  clearly  delineate  the  extent  to  which 
models  based  on  Eqs.  (4)  can  reproduce  the  predictions  of  the 
full  time-dependent  calculation. 

VI .  SUMMARY 

One  conclusion  that  can  be  drawn  from  this  work  is  that 
calculations  of  the  time-dependent  behavior  of  a  junction  using 
the  full  BCS  theory  are  more  feasible  than  one  might  at  first 
suppose,  particularly  if  "memory  reduction"  is  accepted.  This 
is  fortunate  because  as  work  proceeds  on  high-speed  devices 
using  Josephson  junctions,  it  seems  likely  that  it  will  be 
necessary  to  turn  more  and  more  to  this  approach  for  accurate 
predictions.  with  regard  to  memory  reduction,  a  substantial 


* 


-20- 


amount  does  not  change  the  results  by  much,  and  these  changed 
results  are  plausible  and  perhaps  even  closer  to  reality  than 
the  unmodified  case  is. 

As  is  certainly  not  surprising,  the  simple  RSJ  model  is 

shown  to  be  defective  if  quantitative  results  are  desired. 

13 

Another  paper,  dealing  with  single- junction  SQ'JIDs,  shows  this 

even  more  dramatically.  It  deals  with  flux  entry  into  such 

SQUIDs  and  shows  that,  for  some  values  of  the  relevant  parameters, 

the  simple  RSJ  model  is  totally  unsatisfactory. 

Finally,  calculations  such  as  these  can  add  new  insight 

to  our  understanding  of  junction  dynamics.  Equation  (1)  does 

not  lend  itself  readily  to  physical  interpretation,  and  there 

are  undoubtedly  consequences  which  have  not  yet  been  recognized. 

» 

A  minor  example  of  this  is  the  overshoot  that  occurs  in  the 
switching  problem  for  RNC  -  10 r,  and  more  interesting  examples 
may  well  arise  as  work  continues. 
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FIGURE  CAPTIONS 

Figurr  1  Current-voltage  characteristic  for  a  Josephson 

tunnel  junction  with  a  constant  current  source. 

The  line  shows  a  measured  characteristic.  The 

circles  are  calculated  from  the  BCS  theory  with 

R^C  =  lOx.  The  crosses  were  obtained  by  including 

memory  reduction  (see  text)  with  A  *  50.  The 

dashed  line  shows  the  normal  state  result. 

Figure  2  Effect  of  memory  reduction  on  the  coefficients  in 

the  constant  voltage  tunneling  current  equation 

(Eq.  (4)  of  text).  The  solid  lines, taken  from 

ref.  12,  are  the  predictions  of  the  BCS  theory. 

The  points  show  the  effect  of  memory  reduction 

with  the  parameter  A  set  equal  to  50.  Note  the 

rounding  of  the  discontinuities  in  IT-  and  I 

J2  qp 

and  the  fact  that  the  singularity  in  I  ^  is 
replaced  by  a  finite,  broadened  peak. 

Figure  3  Normalized  voltage  versus  normalized  time  for 
switching  to  the  finite  voltage  state  for 
RjjC  =  10t  and  for  a  constant  current  of  1.001  Ic. 

(a)  Predictions  of  the  BCS  theory  and  the  RSJ  model. 

(b)  Prediction  of  the  RSJ  model  with  voltage- 
dependent  resistance.  The  dashed  line  indicates 
the  BCS  result  from  (a)  for  the  longer  times  where 
it  differs.  The  two  are  indistinguishable  on  this 
scale  for  t/x  less  than  about  27,  with  the  oscil¬ 
lations  matching  precisely. 
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Figure  Captions  (continued) 

Figure  4  Switching  for  R^C  =  O.It.  The  RSJ  model  with 

voltage-dependent  resistance/  shown  by  the  dots, 
does  not  agree  as  well  with  the  BCS  prediction  as 
was  the  case  for  R^C  =  10t  (Fig.  3).  However  it 
is  still  much  better  than  the  constant  resistance 
version  (not  shown),  which  does  not  reach  the  first 
peak  until  t/x  =  7.9. 
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Multiple  magnetic  flux  entry  into  superconducting  quantum-interference  devices  (SQUIDs):  A 
general  way  of  examining  the  co*J>  conductance 
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A  new  type  of  experiment  is  proposed  for  obtaining  information  about  the  comb  conductance  of  the 
Josephson  effect  Based  un  measuremem  of  Ruxoid  entry  into  a  superconducting  ring  broken  by  a  Josephsm. 
junction,  the  technique  is  to  operate  in  the  low-damping  regime  for  which  the  voltage  excursions  associated 
with  fluxoid  entry  are  small  For  this  case,  the  constan:  voltage  expression  containing  the  cos4>  conductance 
should  be  valid.  It  is  shown  that  the  erraiicity  associated  with  the  low-damping  regime  has  a  predictable 
statistical  pattern,  which  is  rather  insensitive  to  noise  but  quite  sensitive  to  the  costk  term  A  shunt  resistance 
can  lie  used  to  vary  the  average  voltage  Statistics  .an  be  accumulated  over  a  large  number  of  similar  loops, 
or  over  one  or  a  few  loops  at  slightly  varying  bath  Irnqveranjre  between  runs,  or  even  over  one  lixip  at  one 
temperature  provided  the  noise  Jt  the  junction  has  appropnate  properties  Thus,  the  technique  would  appear 
to  be  capable  ol  estimating  the  controversial  coofUc.ei  i  if  the  cos*  erm  as  a  function  of  voltage  and 
temperature  for  any  type  of  junction  for  which  low  'lamping  can  be  achieved 


I  INTRODUCTION 

The  observation  of  multiple  magnetic  fluxoid 
entry,  or  “quantum  transitions.''  into  superconduct¬ 
ing  quantum-interference  devices  (SQUfDs)— su¬ 
perconducting  rings  which  are  broken  by  Joseph- 
sott  junctions— has  been  reported  in  many  articles, 
usually  with  a  qualitative  or  semiquantitati ve  dis¬ 
cussion  of  the  underlying  reasons  for  the  multiple 
fluxoid  entry  (see,  e.g.,  Ref.  1).  In  a  typical  ex¬ 
periment,  a  magnetic  field  is  applied  in  a  direction 
normal  to  the  plane  of  the  ring.  When  the  applied 
flux  reaches  a  certain  threshold  value,  many  mag¬ 
netic  flux  quanta  will  enter  through  the  junction  in 
rapid  succession.  No  more  fluxoids  enter  the  ring 
until  the  external  flux  reaches  a  second  threshold, 
whereupon  another  group  enters,  etc. 

Some  recent  articles'"  have  reported  results  of 
detailed  simulations  of  this  phenomenon,  using 
simple  models  for  the  junction.  Smith  and  Black¬ 
burn2  have  shown  that  for  rings  with  “high"  damp¬ 
ing  ,iand  Large  values  of  ym  Uc, '+„  (symbols  are 
defined  Liter),  the  number  of  flux  quanta  entering 
the  loop  as  the  external  magnetic  flux  reaches 
threshold  is  a  unique  and  predictable  function  of 
the  SQUID  parameters.  However,  at  small  damp¬ 
ing  the  number  entering  was  shown  by  Wang  and 
Cayley2  to  become  erratic  (or  more  accurately, 
extremely  sensitive  to  small  changes  in  parame¬ 
ters).  The  latter  authors'1  later  inserted  the  cos* 
conductance11'7  into  their  simulations  and  showed 
that  it  can  greatly  affect  the  number  of  fluxoids 
entering  the  loop  in  the  high-damping  region.  The 
measurement  of  the  number  of  fluxoids  entering 
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the  ring  thus  could  give  aseiul  information  about 
the  cos0  conductance,  a  question  which  is  now  very 
much  unsettled.7'1  However,  Cayley  and  Wang1 
cautioned  that  the  high  rapidly  varying  voltages 
expected  to  develop  during  the  flux  entry  might 
render  invalid  the  starting  point  of  the  calcula¬ 
tion,  namely  the  constant-voltage  assumption5'7 
which  leads  to  the  appearance  of  an  explicit  cos* 
term.  That  this  is  probably  true  is  shown  later 
in  the  present  paper,  as  well  as  possibly  by  a  re¬ 
cent  high-damping  experiment5  in  which  the  mea¬ 
sured  flux  entry  did  not  agree  at  all  with  simula¬ 
tions,  with  or  without  the  cos0  term  (see  also  Ref. 

2  in  this  regard). 

In  this  paper  we  consider  especially  the  “errat¬ 
ic”  low-damping  case  in  more  detail,  for  several 
reasons.  First,  erratic  behavior  is  commonly  ob¬ 
served  in  loops’11'2”  and  elsewhere. 21,22  Second, 
as  we  show  later,  the  significant  voltage  chanr-' 
occurs  over  a  time  interval  much  longer  than  the 
inherent  response  times  (picoseconds)  of  super¬ 
conducting  materials,  thus,  validating  the  constant - 
voltage  expression  |Eq.  (1)]  in  the  computations. 
And,  finally,  it  is  evident  that  the  cosf>  term  should 
have  a  marked  influence  on  flux  entry  in  the  low- 
damping  regime  as  well  as  al  high  damping. 

In  the  absence  of  noise  and  for  a  particlular 
choice  of  parameter  values  and  initial  conditions, 
there  is  of  course,  one  value  for  the  final  flux  in 
the  loop.  However,  at  low  damping,  a  small  change 
in  parameter  values  or  a  pulse  of  noise  can  result 
in  a  Lirge  change  in  the  final  flux.  In  this  paper  we 
argue  that  in  spite  of  this,  meaningful  measure¬ 
ments  can  in  fact  be  made  with  such  a  system.  For 
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a  reasonable  amount  of  noise  and  a  reasonable 
range  of  parameter  values,  the  final  flux  will  be 
distributed  in  a  predictable  way  over  a  small  sub¬ 
set  at  values.  Therefore,  the  statistical  distribu- 
tion  at  a  series  of  trials  can  be  predicted.  Since 
the  cos£  term  affects  this  distribution,  the  experi¬ 
ment  should  yield  a  value  for  the  coefficient  of 
this  term.  Further,  since  the  maximum  voltage 
developed  is  controllable  by  a  shunt  resistance, 
such  an  experiment  offers  the  possibility  of  de¬ 
termining  the  cos0  coefficient  as  a  function  of 
voltage,  a  type  of  experiment  not  yet  attempted. 
Tunnel  junctions  are  here  implied,  because  the 
damping  constant  can  readily  be  made  small  with 
them.  However,  the  theory  and  technique  pre¬ 
sented  here  would  apply  to  any  type  of  junction  for 
which  low  damping  could  be  achieved. 

Finally,  in  the  Appendix,  we  show  by  a  simple 
argument  the  source  of  an  empirical  formula  de¬ 
duced  from  the  simulations,2  relating  the  damp¬ 
ing  to  y  for  single  fluxoid  admission.  We  also 
present  a  more  accurate  simple  expression  valid 
to  within  1%  down  to  y  <*  1.  Knowledge  of  this  re¬ 
lationship  is  Important  because  some  Josephson 
devices  based  on  flux  counting  could  give  mis¬ 
leading  results  if  more  than  one  fluxoid  enters  at 
a  time. 

II.  BASICS 

The  circuit  analyzed  is  that  of  the  simple  junc¬ 
tion  shunted  by  capacitance  and  a  phase-dependent 
conductance,  and  connected  to  a  superconducting 
loop  of  inductance  L  (see  Fig.  1).  The  tunnel  cur¬ 
rent  /(/ )  into  the  junction,  indicated  in  Fig.  1,  was 
first  derived  by  Josephson"  for  a  tunnel  junction 
to  be 

Ht)  =  IcBin4>  +  caV+ cos<t>  ,  (1) 

in  which  a  constant  voltage  V  across  the  junction 
was  assumed.  The  coefficients  te,  a0  »  1/R,  and 
a,  are  voltage  dependent.  The  ratio  al/aQ  has  been 
calculated  by  Poulsen15  at  several  temperatures, 
from  tunneling  theory,  and  is  indicated  schemati- 


FIG.  1.  Circuit  used  In  the  simulations  of  this  paper. 
_  The  three  elements  on  the  right  comprise  the  Joseph¬ 
son  junction.  The  tunnel  current  /((),  as  shown  after 
the  capacitance,  is  given  In  Eq.  (1). 


FIG.  2.  Schematic  diagram  illustrating  the  depen¬ 
dence  of  upon  voltage  (solid  curve),  according  to 
tunneling  theory.  The  shaded  rectangle  Indicates  the 
locations  of  experimental  determinations  of  cr i/er  0  for 
tunnel,  microbridge,  and  point  contact  junctions.  A 
recent  experiment  (Ref.  15)  carried  out  on  tunnel  Junc¬ 
tions  at  several  temperatures  near  the  critical  tempera¬ 
ture,  showing  both  signs  for  o,/<7e,  but  opposite  to  that 
expected  from  tunneling  theory.  Is  not  shown. 

cally  in  Fig.  2.  The  rectangle  in  Fig.  2  also  in¬ 
dicates,  with  one  recent  exception,1*  the  results 
of  experiments  designed  to  measure  ajaa,  for 
tunnel,  microbridge,  and  point-contact  junctions. 
The  vertical  dimension  of  the  rectangle  is  meant 
to  suggest  the  error  bars  associated  with  most  of 
the  experiments.  As  is  seen,  the  experiments  tend 
to  agree  with  each  other.  The  results  have  a  sign 
opposite  to  that  predicted  from  tunneling  theory, 
but  are  in  approximate  agreement  with  Landau- 
Ginsburg  theory.1’12  A  recently  published  experi¬ 
ment15  on  tunnel  junctions  at  temperatures  very 
close  to  the  critical  temperature  shows,  however, 
a  sign  change  as  a  function  of  temperature. 

When  an  external  magnetic  field  is  applied  to  the 
loop  of  Fig.  1,  the  differential  equation  describing 
the  circuit  becomes 

+  (1+oi  co8^)  +  2iir  sin^)  +  0  =  ,  (2) 

where  tl  =  t//UC,  P=JLC /RC,  y*Llc/*0,  feisthe 
junction  critical  current  (here  assumed  unaffected 
by  the  magnetic  field  which  would  typically  be 
s  10"*T»  10_aG),  <p  is  the  superconducting  phase 
difference  across  the  junction,  <p,  is  the  applied 
magnetic  flux  in  units  of  ijin,  and  2.068 
x  10~‘*  V s  is  the  flux  quantum.  Ic,  <r0,  and  a,  are 
taken  to  be  constants  in  our  simulations.  The  re¬ 
sistance  R  is  the  quasiparticle  tunneling  resistance 
combined  with  any  shunt  resistance  that  may  be 
present.  In  an  experiment  the  resistance  will  be 
taken  from  the  dc  current-voltage  characteristic 
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of  the  junction.  In  our  calculations,  we  suppose 
that  some  average  R  is  used ,  but  it  would  be  a 
simple  matter  to  use  the  measured  voltage-depen¬ 
dent  resistance  if  this  seemed  desirable.  The  av¬ 
erage  voltage  may  be  appreciably  less  than  the  gap 
voltage,  so  the  tunneling  resistance  may  be  much 
larger  than  the  normal  state  resistance  of  the 
junction.  A  shunt  resistor  may  be  inserted  to  ob¬ 
tain  the  desired  value  of  ^  or  to  adjust  the  value 
of  the  average  voltage. 

Most  of  the  remaining  discussion  in  this  section 
repeats  material  published  earlier. 1,5  We  feel  it 
desirable  to  include  it,  however,  because  of  the 
important  insights  provided  by  the  potential-well 
picture.  Note  that  the  present  0  and  0,  are  equal 
to  the  earlier1,5  ♦  and  ♦,  multiplied  by  2n. 

When  <t>,  is  changing  adiabatically,  Eq.  (2)  can  be 
cast  into  the  intuitively  appealing  form 

ftr.v)=-^(l^cos0)(^)-,  (3, 

where  the  “kinetic  energy"  T  is  ,)2.  and  the 

"potential  energy"  V  is 

V  =  \(0-0,)-2lryCOS0.  (4) 

The  term  ">n  the  right-hand  side  of  Eq.  (3)  repre 
sents  viscous  drag  modulated  by  the  cos<t>  term. 

The  state  of  the  system  may  be  visualized  as  a 
iwrticle  moving  along  the  potential  V  (see  Fig.  3). 
As  the  field  is  raised  slowly  from  zero,  the  par¬ 
ticle  stays  at  the  bottom  of  the  local  well  near  0 
=  0  until  0,  is  large  enough  that  the  local  minimum 
is  now  an  inflection  point  (and  also  the  local  max¬ 
imum),  whereafter  the  particle  begins  to  “slide" 
until  stopped  by  the  damping. 

The  extrema  of  V  are  given  by 

=  0  =  <f>  v  2m  sin0  -  0,  ,  (5) 

■‘0 

and  the  inflection  points  by 

— — i  =0=1  +  2jr>  cos0  .  (6) 

30 

The  solution  of  Eq.  (6)  for  0,  substituted  into  Eq. 


PIG.  a.  "  Potential"  U  nt  fq.  It)  plotted  l<  r  v*a  and 


(5),  specifies  what  0,  must  be  at  the  beginning  of 
the  motion  of  the  particle.  One  finds  for  the 
"break-in”  value  of  the  externally  applied  nor¬ 
malized  flux 

=  2,n +l  t0(>‘2)’  (7) 

and  for  the  starting  inflection  point  <p(, 

0,  -  cos",(-l',2ir>)  =  ’.n  r  l/2rr>  »  0(y~3) .  (8) 

The  next  inflection  point  at  positive  slope  is  of 
course  at  2ir  +  0,.  Its  adjacent  miminum  and  max¬ 
imum  are  at  \nt  (2/> )*/a -*  The  intermedi¬ 

ate  inflection  point  is  at  2ir  -  0f  -  ’  -  \iry  +  0(y~2). 
The  latter  quantities  will  be  used  in  the  Appendix, 
where  wc  consider  threshold  damping  for  single 
fluxoid  admission. 

Figure  3  shows  a  plot  of  the  potential  V,  Eq.  (4), 
for  >  -  5,  and  0,  equal  to  its  “break-in”  value  of 
2n  x  5.2f>3.  Note  tiiai  there  are  about  y  local  wells 
between  the  starting  inflection  point  and  the  bottom 
of  the  "bowl,"  or  overall  potential.  T  hus ,  if  the 
particle  stops  at  the  bottom,  >  fluxoids  have  beer 
admitted  ir.10  the  loop  We  shall  use  the  termi¬ 
nology  “high  damping"  and  “low  damping”  to  mean, 
respectively .  the  cases  in  which  the  particle  slops 
before  reaching  the  bottom,  or  sweeps  past  it. 

Ill  BI-.K  AVION  A I  LOW  DAMPING.  WITHOU I  NOISE 

As  mentioned  earlier,  when  Cayley  and  Wang1 
introduced  the  cos0  term  into  their  simulations  at 
high  damping  and  found  a  large  effect,  they  cau¬ 
tioned  Hint  huge  voltage  excursions  might  be  ex¬ 
pected  which  would  possibly  invalidate  their 
starting  point,  Eq.  (1).  That  this  is  true  is  shown 
in  Fig.  4  for  >  -  100  and  ;?=  12,  a  typical  high-damp- 


HU.  4.  Typical  high-damping  case:  0*12,  V”  100. 
Circuit  parameters  which  give  these  values  of  0  and  V 
a iv,  for  example,  L  *2000  pH,  c’*I4  pl*’,  /t.*0.1  m A, 

H  -  1  Si,  it  [  r  0 .  The  time  and  voltage  scales  are  based 
011  these  v  aloes.  The  gap  voltage  is  not  determined  by 
'tie  /if,,  product,  since  R  is  not  necessarily  the  junction 
Morni.tl  ivmh» 'UK  e . 
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ing  case.  The  voltage  excursions  approach  100% 
of  the  maximum  voltage  developed.  Note  that  there 
is  an  infinitude  of  ways  of  choosing  L,  C,  R,  and 
I c  for  given  >  and  p.  For  the  values  indicated  in 
the  caption,  each  voltage  excursion  takes  place  in 
about  30  ps.  After  fluxoid  entry  begins,  some  <13 
fluxoids  enter  the  loop  in  about  t  ns.  Figure  4 
can  also  be  thought  of  as  particle  velocity  as  a 
function  of  time.  The  particle  is  trapped  0.43  of 
the  way  to  the  bottom  of  the  bowl,  whereupon  it 
undergoes  damped  oscillations  (the  plasma  oscil¬ 
lations)  about  the  bottom  of  the  local  well. 

The  low-damping  regime,  however,  is  quite 
different.  Here  the  voltage  excursions  are  typi¬ 
cally  very  small,  and  only  the  voltage  envelope  is 
important.  Figure  5  illustrates  a  moderate  damp¬ 
ing  case,  with  y  =484  and  fl=  1.77.  Even  on  a 
greatly  expanded  scale  (not  shown)  the  voltage 
excursions  are  scarcely  discernible  until  the  par¬ 
ticle  is  nearly  trapped,  and  even  then  the  ampli¬ 
tudes  -.re  very  small.  For  the  typical  circuit 
parameters  used  in  Fig.  5,  the  voltage  envelope, 
corresponding  to  the  entry  of  about  480  fluxoids, 
develops  over  a  time  of  about  3  ns,  much  longer 
than  the  picosecond  response  times  of  typical 
superconducting  materials.  The  junction  will  thus 
readily  follow  this  adiabatically  varying  voltage, 
and  Eq.  (1)  should  be  a  valid  basis  for  simulations. 
An  accurate  calculation  would  incorporate  the 
voltage  dependence  of  <7,/ou,  but  since  this  is  “un¬ 
known”— this  is  what  the  experiment  is  all  about— 
we  treat  it  as  a  constant;  the  experiment  then 
would  determine  a,/<j0  as  a  function  of  the  average 
voltage  developed  during  flux  entry.  By  means  of 
a  shunt  resistance,  the  latter  can  be  varied  from 
values  comparable  to  the  gap  voltage  to  much 
smaller  values.  In  Fig.  5,  for  example,  suppose 
that  the  gap  voltage  is  2  mV.  Then  the  resistance 
of  2  fl  implies  use  of  a  shunt  resistance  only  slight¬ 
ly  greater  than  2  fl,  since  with  /c  -  0.4  niA  the 
normal  junction  resistance  would  be  about  8  fl. 


FIG.  5.  Moderately  low-damping  case :  0=1.77,  > 

=  4«4.  Circuit  parameters,  which  determine  the  time 
and  voltage  scales,  arc-  L  •  2.r>00  pH,  C-200  pK,  lc 
=  0.4  n'A,  R  =  2  51,  O|»0.  The  dashed  lines  Indicate 
the  voltage  extrema. 


The  resistance  on  the  low-voltage  portion  of  the 
/-  V  curve  actually  being  sampled  would  be  con¬ 
siderably  greater.  As  is  seen,  the  average  voltage 
developed  during  the  flux  entry  is  about  0.2  mV,  or 
10%  of  the  assumed  gap  voltage. 

A  complicating  effect  occurs  in  the  low-damping 
regime,  however.  Wang  and  Gayley1  showed  that 
the  final  state  of  the  system  seems  to  be  erratic. 
Actually,  the  final  state  is  a  very  sensitive  but 
predictable  function  of  the  system  parameters.  An 
example  is  shown  in  Fig.  C  for  p  -  1.20,  o,  -o0,  and 
>  in  the  neighborhood  of  1000.  Here,  there  are  two 
preferred  final  states,  one  near  a  final  flux  num¬ 
ber  of  1100  (the  particle  sweeps  past  the  bottom  of 
the  bowl  and  is  trapped  about  10%  up  on  the  oppo¬ 
site  wall),  and  the  other  near  1000  (the  particle 
does  not  get  trapped  on  the  opposite  wall,  but  falls 
back  to  near  the  bottom).  These  fluxoid  numbers 
are  consistent  with  the  "approximate  theoretical 
maximum  and  minimum"  computed  from  Eq.  (1) 
of  Wang  and  Gayley.3  Gudret23  has  also  found  two 
final  states  in  a  related  calculation.  For  yet  lower 
values  of  damping,  there  can  be  more  than  two 
possible  final  states. 

For  a  damping  constant  near  1.2  and  y  near  1000, 
the  final  fluxoid  number  for  o,  =  +  o0  in,  for  exam¬ 
ple,  the  upper  state  differs  by  only  a  few  from  the 
upper  state  number  for  a,  =  -a0.  Thus,  counting  of 
fluxoids  would  have  to  be  accurate  within  about 
1%  in  this  example,  in  order  to  determine  a  value 
for  ol/oa  with  error  bars  less  than  ±1.  Although 
this  may  be  possible,  there  is  a  much  more  im¬ 
portant  reason  why  fluxoid  counting  from  a  single 
measurement  would  be  inadequate  for  examining 
the  cos0  term  in  low-damj  ng  loops.  As  we  see 
from  Fig.  6,  it  is  not  likely  that  we  could  know  the 
1  of  a  given  loop  to  sufficient  precision  to  predict 
whether  the  final  fluxoid  number  would  be,  for  ex¬ 
ample,  1000  or  1100.  Even  if  we  knew  >  accurate- 
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F  IG.  (i.  Final  llux  values  in  a  superconducting  loop 
with  damping  constant  of  1.20,  n,»a0.  and  y  near  1000. 
Two  states  arc  selected  by  the  system  In  an  almost 
periodic  fashion. 
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ly  on  one  run,  its  value  on  a  subsequent  run  would 
be  slightly  different  if  the  bath  temperature 
changed  slightly.  Thus,  what  would  be  obtained  on 
a  series  of  runs  would  be  a  statistical  distribution 
over  the  possible  final  states. 

We  have  observed  that  the  “occupation  numbers” 
of  the  final  states  are  quite  sensitive  to  the  cos</> 
term,  yet  insensitive  to  noise.  That  is,  the  frac¬ 
tion  of  the  cases  in  which  the  system  will  end  in  a 
given  preferred  final  state  is  se  :sitive  to  the  value 
of  c,/o0,  and  this  fraction  is  readily  determined 
from  the  simulations  by  varying  >  in  small  incre¬ 
ments.  Thus,  instead  of  having  to  count  fluxoids 
with  great  precision,  one  has  the  attracts  ~  al¬ 
ternative  of  accumulating  statistics,  using  rela¬ 
tively  crude  fluxoid  counting.  In  such  a  proce¬ 
dure,  one  could  fabricate  a  series  of  supercon¬ 
ducting  rings  with  closely  similar  values  of  y. 

The  individual  values  would  presumably  be  dis¬ 
tributed  somewhat  randomly  witliin  a  small  range, 
and  the  results  of  the  flux  measurements  would  be 
a  statistical  distribution  of  values  among  the  pre¬ 
ferred  final  states. 

Rather  than  making  measurements  on  a  Large 
number  of  similar  loops,  one  would  probably  pre¬ 
fer  to  make  repeated  measurements  on  one  or  a 
few  loops,  deliberately  varying  the  bath  tempera¬ 
ture  slightly  between  measurements,  for  example, 
a  few  millidegrees.  Since  y  is  proportional  to  le, 
which  in  turn  is  proportional  to  the  temperature- 
dependent  energy  gap,  this  procedure  would  seem 
to  be  a  convenient  method  for  varying  >  in  small 
steps.  Still  another  procedure  (here  anticipating 
the  results  of  Sec.  IV)  might  be  to  make  repeated 
measurements  on  one  or  a  few  junctions  at  a  “fixed” 
temperature,  allowing  the  noise  at  the  junction  to 
be  the  statistical  generator.  As  we  shall  see  in 
Sec.  IV,  the  noise  parameters  would  have  to  fall 
in  a  certain  range. 

Table  I  shows  the  results  of  a  ratio  analysis  at 
,)  =  1.18,  1.20,  and  1.22,  which  is  a  reasonable 
range  of  uncertainty  for  this  parameter.  The 
Table  entries  show  the  percentage  of  the  total  num¬ 
ber  of  possible  times  that  the  flux  entry  will  be 
found  in  the  upper  of  the  two  final  states;  they  are 
determined  by  incrementing  >  in  very  small  steps, 
and  are  accurate  to  *  0.1  (for  the  noise-free  en¬ 
tries).  The  differenre  of  more  than  a  factor  of  2 
between  the  results  at  ct,  ct0  -  ♦  1  and  -1  should 
make  the  distinction  between  these  values  readily 
discernible.  The  ultimate  accuracy  of  the  deter¬ 
mination  of  o,  o.,  will  depend  principally  upon  the 
amount  of  statistics  accumulated.  The  effects  of 
noise  are  discussed  in  Sec.  IV. 

There  are  of  course  other  values  of  Hand  >  which 
may  be  appropriate.  We  selected  the  range  dis¬ 
played  in  Table  I  because  just  two  well-separated 
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TABLE  1.  The  entries  in  the  e,  columns  are  the  per¬ 
centage  of  occurrences  for  which  the  final  flux  entry  into 
the  superconducting  loop  is  in  the  high-flux  state  (about 
1100  in  this  case).  Values  of  y  near  1000  are  used,  with 
damping  values  fi  as  shown.  The  last  row  indicates  the 
results  of  noise  simulations  at  a  noise  frequency  of  1000/ 
./.7s  and  a  noise  amplitude  of  O.lrt,.  Some  228  simula¬ 
tions  using  different  sets  of  (pseudo)  ramdom  numbers 
are  used  lor  each  value  of  n,  in  the  noise  simulations. 

The  uneertaintieo  indicated  correspond  to  one  standard 
deviation  (lief.  2j>. 
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final  slates  occur.  Lower  damping  might  result 
in,  for  example,  three  final  states,  two  of  which 
would  likely  be  fairly  close  to  each  other,  and  pos¬ 
sibly  confuse  the  results.  However,  this  is  not 
necessarily  so,  because  one  might  simply  examine 
the  well-separated  state,  unless  its  percentage 
occupancy  is  so  low  that  adequate  accumulation  of 
statistics  becomes  difficult. 

IV  EFKECrSOI  NOISE 

It  is  important  to  understand  the  effects  of  noise 
on  the  results  of  Sec.  111.  Noise  may  be  pictured 
as  a  rocking  of  the  "bowl,"  representing  the  po¬ 
tential  energy  V,  about  its  bottom.  That  is,  from 
Eq.  (4),  the  variation  of  the  potential  due  to  a 
variation  A0,  in  <px  is 

&V  -  [<px  -  ip)(i<px .  l9) 

Thus,  one  might  anticipate  that  noise  would  tend 
to  make  the  flux  entry  more  regular— the  particle 
tends  to  get  shaken  toward  the  bottom  of  the  bowl. 
For  a  particle  trapped  high  on  the  wall  of  the  bowl  . 
which,  can  occur  only  for  very  low  damping,  noise 
cun  fairly  easily  displace  the  particle  from  its 
relatively  shallow  local  well.  Bu*  for  moderately 
low  damping,  in  which  the  particle  is  trapped,  for 
example,  10',',  above  t lie  bottom  of  the  bowl,  noise 
of  “ordinary"  amplitude  will  not  displace  the  par¬ 
ticle  once  the  particle  lias  settled  down,  which 
occurs  after  several  plasma  oscillation  periods. 
One  should  also  note  that  in  addition  to  tending  to 
prevent  trapping  in  a  certain  well,  noise  can  also 
promote  trapping  in  a  well  in  which  the  particle 
would  not  remain  in  the  noise-free  case.  There  is 
a  strong  symmetry  in  these  two  cases,  which  our 
simulations  bear  out. 

Noise  enters  primarily  through  o,—  even  Johnson 
noise  in  the  junction  or  shunt  resistance'"  can  be 
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lumped  into  <pr.  A  “second-order”  effect  would  lie 
in  the  variations  of  >  because  of  variations  in  4. 
e.g.,  because  0,  is  noisy,  but  this  can  surely  be  ig¬ 
nored.  Thermal  fluctuations  could  also  affect  4 
through  the  energy  gap,  but  both  the  amplitude  and 
frequency  would  be  too  low  to  be  of  significance 
for  temperatures  not  too  close  to  Tc. 

We  have  examined  the  effects  of  noise  by  adding 
to  0,  a  suitably  distributed  (pseudo)random  num¬ 
ber  at  each  new  time  increment  in  the  computation 
(or  at  some  multiple  thereof,  thereby  varying  the 
noise  frequency,  or  spectrum).  Noise  levels  (stan¬ 
dard  donations)  up  to  1CK?  of  0,  have  been  used. 

It  is  evident  that  the  “erraticity"  at  low  damping 
will  be  dependent  upon  noise  frequency.  If  the 
particle  has  not  had  time  to  settle  to  the  bottom  of 
a  local  well,  a  pulse  of  noise  has  a  fair  probability 
of  kicking  the  particle  out.  Thus,  the  natural  fre¬ 
quency  against  which  to  compare  noise  frequency 
is  the  plasma  frequency,  which  is  the  frequency  of 
the  motion  of  the  particle  in  a  local  well.  Ttiis  is 
readily  obtained  from  Eq.  (2)  by  setting  <p  -  Inn 
♦  <S>'  where  >0'  I  c<  1  and  n  is  an  integer.  The  ho¬ 
mogeneous  solution  to  the  equation  thus  linearized 
(and  here  dropping  the  “cos 0  term")  is  exp (/>/,) 
where 

P=-\„i\ii2~H\*2ny)Y'\  (10) 

Since  ft1  is  less  than  about  9>  for  multiple  flux 
entry  (see  the  Appendix),  this  is  an  underdamped 
case.  Further,  if  /32  <;  8 iry  ,  which  is  the  ordinary 
case  for  low  damping,  the  angular  plasma  frequen¬ 
cy  w,  is  just  (8iry)1/s  in  1 1  space,  or  (81 ry/ LC)'‘n  in 
real  time.  The  damping  time  is  2/;je  jy‘/s,  the 
inequality  being  the  condition  for  multiple  transi¬ 
tions,  as  above.  Thus,  the  damping  time  is  great¬ 
er  than  the  plasma  period  2/(8>  'ir)'n,  and  the  par¬ 
ticle  will  always  make  several  swings  through  the 
bottom  of  a  local  well  before  coming  to  rest. 

The  plasma  period  has  a  significance  beyond 
that  of  a  trapped  particle:  a  particle  traveling 
slowly,  but  not  quite  slowly  enough  to  be  trapped, 
clearly  will  require  about  one-half  the  plasma 
lieriod  in  traveling  from  one  local  maximum  to  the 
next.  Thus,  we  can  easily  estimate  the  time  re¬ 
quired,  or  the  number  of  noise  pulses  oecuring, 
when  the  particle  moves  from  a  fluxoid  number 
of,  for  example,  1100  to  1000,  to  use  the  example 
of  Sec.  ID. 

There  are  two  effects  ol  noise  which  must  be 
considered,  u)  By  tv  w  much  does  noise  produce 
scatter  about  the  preferred  final  states,  and  arc 
the  mean  values  shifted  from  the  noise-free  ease'' 
oi)  By  how  much  does  noise  affect  the  ration  of 
the  “ocrupation  numbers”  in  the  (now  spread-out) 
final  states?  We  anticipate  that  the  preferred  ex¬ 
perimental  procedure  would  be  that  of  determining 
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FIG.  7.  Histograms  showing  how  different  non  e 
frequencies  affect  the  final  flux  values  In  a  superconduct¬ 
ing  loop.  Noise  amplitude -O.u.'i  if  =1.20;  u  ,  =  0, 
nine  values  of  y  in  the  range  from  999.0  to  Iflol.7. 

Alv.it  SO  runs  a  ere  made  at  each  noise  frequency. 


(lie  ratios.  Thus,  question  (i)  is  not  important  un¬ 
less  the  scatter  is  so  great  that  the  clusters  over¬ 
lap. 

That  t he  mean  values  of  the  final  states  should 
not  be  materially  affected  by  noise,  as  seen  in 
Fig.  7,  is  intuitive:  The  preferred  states  are  de¬ 
termined  by  the  values  of  p  and  > .  If,  for  example, 
a  small  change  of  >  changes  the  trapping  point 
from  the  upper  to  the  lower  of  two  preferred 
states,  by  the  same  reasoning  a  noise  kick  on  the 
trapped  particle  near  the  extreme  of  its  plasma 
oscillation  swing  can  e)ect  the  particle  and  cause 
it  to  seek  the  lower  preferred  state. 

Scattering  obviously  increases  with  increasing 
noise  amplitude.  The  example  of  Fig.  7  shows  no 
overlap  between  clusters  even  at  the  relatively 
large  noise  amplitude  of  0.050,.  The  influence  of 
noise  frequency  upon  the  scatter  is,  however,  not 
so  easy  to  understand,  and  certainly  not  easy  to  ex¬ 
plain  in  a  few  words.  Since  the  question  is  not  impor¬ 
tant  m  the  present  context  we  shall offc-  only  the  fol¬ 
lowing:  When  the  noise  frequency  is  low,  Lhelast 
noise  pulse  before  trapping  (or  nontrapping)  is  im¬ 
portant  and  may  cause  trapping  in  an  "abnormal" 
state,  thus,  causing  scatter  When  the  noise  fr<  - 
quency  is  high,  however,  there  is  a  great  amount  of 
sell  cancellation  of  the  noise,  and  the  scatter  is 
smaller.  Our  s  imu  la  turns  do  show  that  for  noise  fre¬ 
quencies  w,  less  than  ur,  there  is  a  good  tut  ol 
scatter  about  the  final  states,  when  as  for 

the  clustering  is  quite  narrow.  Figure  7  il¬ 
lustrates  tlus  fir  three  values  of  noise  frequency; 
10  ■,  TT~ ,  100  1  ,  and  1000  ,T~  .  These  are  to 

be  compared  to  the  plasma  frequency  ,7  7  -  tV 
'  25  .  X  for  >  1000. 

The  important  question  is  number  Uil  above. 

One  expects  that  the  occupation  number  ratios 
will  not  be  greatly  affected  by  "reasonable"  noise 
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amplitudes  since  the  trapping  and  nontrapping  sit¬ 
uations  should  be  affected  rather  symmetrically 
by  the  noise.  All  the  simulations  we  have  made 
bear  this  out.  Table  I  illustrates  this  for  >  near 
1000  and  ,3  near  1.2,  at  the  very  large  noise  am¬ 
plitude  of  0.1  4>x  and  a  noise  frequency  of 
1000,  iiC.  The  228  calculations  made  for  each 
value  of  Oj  result  in  an  uncertainty  in  the  listed 
percentages  as  shown  by  the  indicated  standard 
deviations. 25  Calculations  at  smaller  noise  am¬ 
plitudes  have  also  been  made,  and  also  show  the 
insensitivity  of  the  percentages  to  noise. 

Finally,  our  simulations  with  noise  show  tliat  the 
percentages  are  not  biased  by  choice  of  } ,  for  a 
given  J.  For  example,  if  we  consider  only  those 
values  of  > ,  within  the  range  999-1001,  for  which 
the  final  state  would  be  the  lower  of  the  two  pos- 
siblestates  in  thenoise-free  ease,  we  find  the  same 
percentages  as  given  in  Table  I  to  within  our  sta  ¬ 
tistical  significance.  Thus,  an  experiment  making 
repeated  runs  on  a  single  junction,  allowing  the 
noise  at  the  junction  to  generate  the  statistics, 
may  be  an  acceptable  technique.  Of  course,  the 
noise  amplitude  has  to  be  large  enough  to  be  ef¬ 
fective,  but  not  so  large  as  to  produce  chaotic  re¬ 
sults.  Our  simulations  have  shown  that  quite  a 
braid  range  of  noise  amplitude  is  acceptable.  A 
noise  frequency  comparable  to  or  larger  than  the 
plasma  frequency  would  also  be  desirable. 


V.  CONCLUSIONS 

We  conclude  that  it  should  be  possible  to  derive 
useful  information  about  the  cos0  conductance  by 
measurement  of  fluxoid  entry  into  a  supercon¬ 
ducting  loop  at  low  damping.  Low  damping  ensures 
that  the  basic  equation  containing  the  cos0  term  is 
a  valid  basis  for  this  type  of  experiment. 

The  procedure  described  in  this  paper  should 
greatly  improve  our  knowledge  of  a.  at  low  vol¬ 
tage.  Moreover,  by  use  of  a  shunt  resistance,  the 
average  voltage  developed  during  flux  entry  can 
be  controlled,  and  thus  a,  a„  can  be  estimated  as 
a  function  of  voltage,  a  dependence  which  lias  not 
been  measured  to  date.  Since  tunneling  theory 
predicts  a  large  discontinuity  in  e,  <j,,  at  the  gap 
voltage,  with  a  sign  reversal,  measurements  near 
such  voltage  would  be  particularly  exciting. 

The  temperature  dependence  of  o,  o,,  could  also 
bo  obtained,  obviously,  by  varying  the  bath  tem¬ 
perature.  It  would  be  interesting  to  compare  such 
results  with  those  recently  obtained.15 

Noise  is  not  likely  to  be  an  obscuring  factor, 
according  to  our  calculations  and  may  be  desirable 
in  accumulating  statistics. 


Finally,  provided  low  damping  cun  be  achieved, 
the  experiment  can  be  performed  with  any  type  of 
Josephson  junction.  This  could  be  very  interesting 
because  to  date  each  type  of  junction— tunnel,  mi¬ 
crobridge,  or  point  contact— has  used  an  entirely 
different  experimental  method  for  examining  the 
cos0  conductance. 

Measurement  of  fluxoid  entry  at  hi^h  damping 
can  also  be  made,  but  as  Cayley  and  Wang1  pointed 
out,  a  theory  allowing  for  a  dynamic  voltage28 
should  then  be  used  to  calculate  the  expected  flux¬ 
oid  number.  Since  the  theory  for  large  rapidly 
varying  voltage  would  not  have  an  explicit  cos0 
term,  one  could  not  then  speak  of  determining  the 
value  of  j  „  bat  only  of  confirming  or  denying 
tin-  validity  of  the  complete  theory. 


A(k  NOW  U.DGME.NTS 

We  wish  to  thank  our  colleagues  in  the  cryoelec¬ 
tronics  group  of  the  Natl.  Bur.  Stand.  (U.  S.)  for 
many  invigorating  discussions.  Part  of  this  work 
was  done  while  one  of  us  (R.  I.  G.)  was  a  guest 
worker  at  Natl.  Bur.  Stand.  (U.  S.).  He  wishes  to 
thank  the  staff  for  its  hospitatily.  This  work  was 
supported  in  part  by  the  Office  of  Naval  Research 
Contract  No.  N0001477C0415. 


APPENDIX 

In  tlus  Appendix  we  consider  the  threshold 
damping  for  single  fluxoid  entry,  in  order  to  show 
the  source  of  a  relation  deduced  empirically  from 
earlier  simulations. We  also  develop  a  more 
accurate  expression  valid  over  a  larger  range  of 

To  examine  the  threshold  damping  in  question, 
we  ask:  What  must  be  the  value  of  the  damping 
constant  .1  in  order  tlut  the  particle  not  slide  be- 
yo.-.d  the  first  maximum  ol  the  potential  TJ  at 
s  •  (2  >)'  :?  (The  positions  of  the  first  lew  ex¬ 
trema  and  inflection  points  of  ’ll  were  obtained  in 
Sec.  1J.)  The  rotation  deduced  empirically'’5  for 
;  in  the  region  above  100  is  ,i,  -3,0>‘'2.  To  show 
the  source  of  this  relation  we  use  the  following 
simple  argument.  First,  we  observe  that  for  7 
»>  1.  as  illustrated  in  Fig.  8  for  7  =200,  the  poten¬ 
tial  energy  curve  is  quite  steplike  initially,  rather 
than  having  pronounced  minima  and  maxima.  This 
suggests  approximating  the  potential  by  the  piece¬ 
wise  linear  portions  shown  in  Fig.  8.  We  draw  a 
horizontal  line  through  0  ,..  and  a  straight  lino 
with  the  correct  negative  slope  of  -4ri  through  ttie 
inflection  point  near  <;i  '.r.  One  can  easily  show. 
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FIG.  8.  First  portion  of  the  potential X)  for  y  =  200  ami 
<h,  =  2*  >2002250,  together  with  the  piecewise  linear  portions 
used  to  approximated  in  this  region.  The  vertical  scale 
is  arbitrary. 


once  fa  is  estimated,  that  the  particle  sliding  down 
the  linear  slope  will  nearly  have  its  terminal  ve¬ 
locity 

dp 

<"  i 


(Al) 


when  the  horizontal  portion  is  reached.  On  a  hor¬ 
izontal  line,  the  particle  travels  a  distance  equal  to 
its  initial  velocity  divided  by  the  damping  constant; 
that  is, 


0* 


1  dtp 

'  ,3  til , 


(A2) 


Setting  the  asymptote  <pKor  equal  to  §»,  and  tp^,, 
equal  to  2ir  [the  intersection  of  the  two  straight 
lines  is  at  2#  +  CKy-*)],  requires  fa  =  (2/ir) 
x(dtp  dt  i)lltimM„  which  with  Eq.  (Al)  gives 

A=(8y)‘/2,  (A3) 

which  has  the  observed  square-root  behavior. 

Note  that  8l/I  =  2.83.  Note  also  that  if  in  Eq.  (A2) 
the  asymptote  is  refined  to  ’ir «  (2/y)l/*  then  Eq. 
(A3)  becomes 

*-(M‘/*-4/»;  (A4) 

that  is,  an  offset  is  predicted.  We  have  made  a 
computer  study  of  <3,  vs  >  down  to  values  for  which 
jti,  vanishes.  The  expression 

2.99yl/2  -  2.53  ,  (A5) 

is  accurate  to  within  1%  for  all  values  of  >  above 
unity.  For  y  <1,  fa  dips  slightly  below  this  line 
and  falls  to  zero  at  y,ml,  =0.733.  At  this  value,  the 
second  maximum  of  v  has  risen  to  the  level  of  the 
starting  point,  so  that  even  at  zero  damping,  not 
more  than  one  fluxoid  can  enter  the  loop.  is 

the  solution  of 

2ir>  eos(4ir2>2  -  l)1'2  =  -1  ,  (A6) 

which  results  from  requiring  that  tp  *  2 try  sin0 
have  the  same  value  at  0,  as  at  4a -0,. 

Knowledge  of  threshold  damping  is  important  be¬ 
cause  any  Josephson  device  acting  as  a  flux  coun¬ 
ter  in  some  measurement  process  could  give  quite 
misleading  results  if  more  than  one  fluxoid  would 
enter  at  a  time.  For  example,  in  the  analog-to- 
digital  conversion  of  a  continuous  signal  which 
causes  fluxoid  entry  or  expulsion  from  a  super¬ 
conducting  loop,  the  reconstructed  signal  could 
be  quite  distorted  if  this  occurred.27 
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